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I. INTRODUCTION 



Parton distribution functions (PDF)s contain nonperturbative information about the in- 
trinsic structure of hadrons in terms of their constituents — quarks and gluons [l|, Q, [sj. 
In completely inclusive processes (e.g., deeply inelastic scattering (DIS), where the hard 
virtual photon with momentum = —Q^ probes a hadron h with momentum P), inte- 
grated PDFs fi{x, Q"^) {i marking the sort of parton) depend on the longitudinal- momentum 
fraction x, which becomes equal to the Bjorken variable xb in the limit — > oo, 
Xb = Q^/2(q' ■ P) = const, and on the scale of the hard subprocess Q^. In the Bjorken 
limit, these distributions (Feynman parton densities) are related to the (unpolarized) quark 
and antiquark structure functions 



(1) 



in the leading-twist approximation and in leading order of the coupling (where is the 
electric charge of the quark of flavor i) j^]. Equation (p[]) originates from the DIS factorization 
expression 



F,{XB,Q') = J2 f- 



Xb Q"^ . . 

z 



(2) 



where the perturbative coefficient functions Cu are taken in leading order (LO): Cif 



1^2 



ef 5{x/z — 1). It can be considered as a triumph of perturbative QCD that the QCD 
evolution correctly describes the logarithmic dependence of the parton densities on the hard 
scale and is, therefore, able to explain the experimentally observed violation of the Bjorken 
scaling. Moreover, the gauge-invariant operator definition of integrated PDFs 



fi/h{x,Q'^ 
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-ig f dz^^A^izY'^ 




Jo- 





(3) 



with A;+ = xP^ allows one to relate their moments to the matrix elements of the twist-two 
operators arising in the operator product expansion (OPE) on the light-cone J5|. The renor- 
malization properties of these PDFs are governed by the DGLAP equation 0, [gJ , establishing 
the logarithmic dependence on mentioned above. 

The study of semi-inclusive processes, such as semi-inclusive deeply inelastic scattering 
(SIDIS), or the Drell-Yan (DY) process — where one more final or initial hadron is detected 
and its transverse momentum is observed — requires the introduction of more complicated 
quantities, viz., unintegrated, i.e., transverse-momentum dependent (TMD) distribution or 
fragmentation functions: 



semi-inclusive 



fi/h{x, k^,Q^ 



(4) 



The most natural generalization of the operator definition (j3]), respecting gauge invariance 
and collinear factorization, reads 0, 0, Isl, 0, [13] 



fi/h {x, kj_; 11^) 



2 J 27r(27r)2 

X [oo~, oo", oo J^7+[oo", oo^; oo~, J [oo~, 0^; 0", J^i(0", 0±)\h) 



(5) 
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where gauge invariance is ensured by means of the path-ordered Wilson-hne operator (gauge 
hnk) with the generic form 



[y, x\V] = V exp 



x\T] 



(6) 



The transverse gauge hnks extending to hght-cone infinity are also included in ^ and the 
dependence on is taken into account via the renormalization-group (RG). 

However, as it has been pointed out in [3], when one retains in the parton densities the 
intrinsic transverse momentum, extra undesirable divergences appear. These divergences are 
associated with the particular features of the light-cone gauge (or the use of purel y li g ht like 
Wilson lines) that must be removed by some consistent method (see, e.g., jl, El, 12, 13, 14]). 
For instance, they can be avoided by using non-lightlike gauge links in covariant gauges, or by 



employing an axial gauge, but going off-the-light-cone [7|, ll5||. This involves the introduction 



of an extra rapidity parameter and entails an additional evolution equation [3], rendering 
the reduction to the integrated PDF questionable. 

Another strategy, based on a subtraction formalism of these extra divergences in terms 
of a "soft" factor (defined as the vacuum average of particular Wilson lines and amounting 
to a generalized renormalization of TMD PDFs), was presented in Refs. (l3l . 



The major finding in our previous investigations in [l3|, 3] was that, adopting the light- 
cone gauge, the leading gluon radiative corrections associated with the transverse gauge link 
were found to give rise to an extra term in the anomalous dimension of the TMD PDF that 
exhibits a lnp+ behavior — characteristic of a contour with a cusp. The new definition for 
the TMD PDF, proposed in these works, has two important advantages: (i) it reduces to the 
correct integrated case and (ii) it coincides with the result obtained in the Feynman gauge 
that is untainted by contour obstructions. 

In (non-covariant) axial gauges the partonic interpretation of the distribution functions is 
preserved because the gauge links can be set to unity under the gauge condition. It is well- 
known that among the axial light-cone gauges there is one which has important advantages. 
Indeed, employing the light-cone gauge in association with the Mandelstam-Leibbrandt (ML) 
prescriptio n |22l . |23|^ in the calculation of the quark self-energy and the quark-quark-gluon 
vertex fiol . |20|. and also the DGLAP kernel in NLO [2l[, it was shown that no undesirable 
singularities appear — even at the intermediate steps of the calculation — in contrast to the 



principal- value (PV) prescription used in [2J]. In the ML-gauge, the contributions of the real 
and the virtual diagrams are well-defined in the "end-point" region separately. [This is the 
region proportional to the delta-function of the longitudinal fraction of the hadron's momen- 
tum ~ 5(1 — x)]. In the calculation of the evolution of the inclusive (integrated) PDF, this 
is not important — at least in LO — since all these singularities cancel in the final result. But 
the situation changes for the unintegrated TMD PDFs. In that case, the spacelike distance 
between the quark operators in the corresponding matrix elements acts like an ultravio- 
let (UV)-regulator, thus preventing the mutual cancelation of the extra (mixed) ^ divergent 
terms between the virtual and the real gluon contributions. As a result, the rapidity diver- 
gences contributing to the UV-divergent part of the self-energy graph remain uncanceled. 
Exactly those terms of the splitting function, arising from the "end-point" region, give rise 



^ Let us call in what follows the light-cone gauge with this prescription the "ML-gauge" . 
^ We distinguish between "pure" singular terms, having only a single UV or rapidity divergence, and "mixed' 
ones, which contain two poles of different origin simultaneously. 
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FIG. 1: Integration contour and poles in the (Re q^, Im q^) plane: the poles of the gluon propagator 
using the ML-prescription (position 1) and those in a covariant gauge (position 2) belong to the 
same, i.e., second and fourth, quadrants. This is in contrast to the poles pertaining to the principal- 
value prescription (position 3). The Wick rotation can be performed without changing the position 
of the poles. 



to extra (mixed) UV- singularities in the TMD function that can be eliminated by employing 
the ML-prescription — even extending the calculation to the NLO 211 . 



From the field-theoretical point of view, the ML-gauge has very attractive properties as 
well. Due to the position of the poles in the same quadrants as for the free gluon propagator 
(see Fig. [1]), one can readily perform the Wick rotation to the Euclidean space. This is not 
possible for the ^"-independent prescriptions. Thus, in the ML-gauge, the standard power 
counting rules allow one to estimate the UV-divergences. Moreover, is has been shown that 
the ML-prescription arises naturally 

We have pointed out recently [13 



14l | that the renormalization-group properties of TMD 
PDFs can be profitably analyzed in terms of their UV anomalous dimensions. The main 



in a consistent quantization procedure [25|, l26|, [27 



reason is that anomalous dimensions are local quantities originating from the geometrical 
obstructions of the gauge contours: endpoints, cusps, or self-intersections. Within this con- 
text, gauge-invariant quantities have to fulfil anomalous- dimension sum rules that represent 



logarithmic, i.e., additive, versions of the Slavnov- Taylor identities [IJ]. In our previous 
works all the calculations have been done in the light-cone axial gauge with additional q~- 
independent pole prescriptions, notably, the advanced, retarded, and the principal-value one. 
It was shown that the extra contribution to the anomalous dimension of the TMD PDF, 
given by Eq. ([2]), can be identified (at least in the one-loop order) with the well-known 



cusp anomalous dimension [28!]. In the present investigation we apply this type of approach 
to a pole configuration of the gluon propagator controlled by the ML prescription. In what 
follows, we shall first analyze (Sec.[Tll) the behavior of the transverse component of the ("clas- 
sical") gauge field at light-cone infinity — required for the derivation of the transverse gauge 
link that eliminates the residual gauge freedom (after fixing the gauge by A'^ = 0) — and 
derive an explicit expression for the gauge field in the ML-gauge. Then, we shall calcu- 
late the UV-divergent parts and the corresponding anomalous dimension of the TMD PDF 
(Sec. imp . It is remarkable that the result obtained this way is free of undesirable terms 
related to contour obstructions and coincides with the double anomalous dimension of the 
fermion field. As we shall show later (Sec. lIVp . the crucial point in verifying the validity 
of the modified definition of TMD PDFs, Eq. (1271) below, proposed in 13|, [ij] , is that the 
so-called mixed rapidity divergences are absent in the ML-gauge, while the contribution of 
the soft factor (which has been introduced in order to cancel these divergences in the case 
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of q -independent prescriptions) reduces to unity in the one-loop order — see Sec. |Vl Our 
conclusions are presented in Sec. |Vll 



II. TRANSVERSE GAUGE FIELD AT LIGHT-CONE INFINITY IN THE ML- 
GAUGE 

It was argued in Refs. [sl, [9| that in order to restore the contribution of the gluon 
exchanges between the struck quark and the spectator in the light-cone gauge, one has to 
take into account the accumulation of the corresponding phase in the transverse direction 
[Note that this phase is suppressed in covariant gauges] . It is precisely this phase that yields 
the additional transverse gauge link at light-cone infinity, introduced in Refs. [sl, 0, [lO] , the 
reason being that this phase is accumulated very slowly. Below, we derive an expression 
for the gauge field in this situation by adopting the ML-prescription. This has not been 
considered before in the literature. 

We commence our analysis by calculating the gauge field, the source of which is a "clas- 
sical" current 

3i.{y) = 9 j dy'^ 5^^\y - y') , y'^ = v^T , (7) 

corresponding to a charged point-like particle (e.g., a struck quark in a SIDIS process) and 
moving with the quasi-constant four-velocity along the straight line w^r. Note that the 
velocity changes only at the origin, where the sudden collision with the hard photon takes 
place and the quark is derailed to its new "trajectory". The gauge field related to such a 
current is given by 

A^{i)= !d'yV^''{i-y)3M , (8) 



where D^^'^ is the gluon Green's function. We assume that the velocity of the struck quark 
is parallel to the "plus"- and the "minus"- light-cone vectors before and after the hard 
collision, respectively: 



n 
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n^r) + n, 



dr 6^^\ 



dq 



-iq y 



27r q- + iO 



- nj{y^) 



y — n T 
dq~^ e 
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iO 



Using the results of Refs. |l3l . Il4l | , we write 

dS 



-9 



dy+dy-d'y^e'^-y6{y~)6^'\y^) 



2(27r)4 

where the free gluon propagator in the light-cone gauge = has the form 



d\ 
(2^ 



-iq-z 



dS _ 
(2^ q^ 



-iq-z 



qf'in-y + q^in-y 



(9) 



(10) 



with the square bracket being used in order to remind that [q'^] has yet to be defined. Here, 
and in what follows, we neglect the quark and gluon masses m and A, since we are mainly 
interested in the UV-singularities. 
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Performing the integration over the variable we get 

^.(,) . / J?X e-''< 4a (r - I . (12) 



2(27r)4 q^ + iO V [g 

Before continuing, let us first verify that the longitudinal (light-cone) components of the 
gauge field vanish and that there is no contradiction with the gauge condition. The "plus"- 
component reads 

■^^-(^-'v-)-;(("*r-^^if^P)-i-|f:j. (13) 

Notice that one has for both types of pole prescriptions: q~-independent, as well as for 
^"-dependent ones (like the ML-prescription), q~^/[q~^] = 1 (see, e.g., Ref. [29|), so that 

A+ = , 

a result in agreement with the light-cone gauge. One the other hand, for the "minus" - 
component, we start with 

■4- - ■ <) - < (("^)'' - ^^^f^) - - ^ . (") 

and carrying out the integration over q~ in Eq. fll2p with 5{q~), we find 

yl^ = . 

Now turn to the evaluation of the transverse components. Let us recall that employing 
^"-independent prescriptions, the transverse component of the gauge field is 

A^(oo-^^) = |-Coo V^ln Al^^l, (15) 

where the numerical constant Coo depends on the pole prescription according to (see [9| 

, Advanced 

Coo= { -1 , Retarded (16) 
-\ , Principal Value . 

Here A is an IR-regulator that does not enter the final results. 

The ML-prescription, being dependent on both variables g"*" and g~, gives rise to a more 
complicated pole structure in the complex g° plane, viz., 

1 



ML 



q++iOq 1^27^ 



q+q~+iQ 



The two possible forms of this prescription, displayed in Eq. f[T7|) . are, in fact, equivalent to 
each other. After performing the integral over y in Eq. ( |T0|) . one finds 



2(27r)4 q^ + iO y ' \q 



ML 
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Taking into account that the hght-cone vectors have only longitudinal components and 
separating out the transverse integrations, we get 



- « ^ y (# «^ e-^ J ^ .(,-) (^.-5)^ ^ (19) 

In order to compute the longitudinal part of this expression, we use for the denominator 
the a-representation and employ for the delta-function the integral representation 

This allows us to write 

dq^dq- e-iii+r+'i-i+) fdq+dq- e-ii<i^ 



(27r)2 ^ (g2 + zO)[g+]ML J (27r)2 ' (q^ + tO)[q+]ML 



(2vr)3 7o y [g+]ML 

X exp [i{aq^ + (^+ + A)g~ + T?^)] • (20) 



To proceed, we make use of the results obtained in (that can be directly derived 

by applying Cauchy's theorem) 



-I 

a 



(21) 



f dq^dq r-/ + - n+ ~ n~ 4-M 27r f / 

and employ the representation 

^ f 1 - e-'«"^) = %C [ dr e"'"«"^ (22) 



to find 

Evaluation of the pure transverse integral yields 

where A is again an auxiliary IR regulator which shall ultimately drop out from all physical 
quantities. Therefore, the transverse gauge field at light-cone infinity in the ML-gauge reads 

A^(oo-;^J = -^V^ln A|^^| . (25) 
47r 

Hence, imposing to the evaluation of the gluon propagator the ML-prescription, the trans- 
verse gauge field at light-cone infinity is given by Eq. ( l25l) and is a total transverse derivative, 
just as its counterpart f|T5l) in the case of ^"-independent prescriptions. But there is a crucial 
difference: in contrast to a ^"-independent prescription, the ML result does not bear any 
dependence on the imposed boundary conditions encoded in the constant Coo- The result- 
ing phase, accumulated by the struck quark moving along the "plus" -light-cone ray, can. 
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therefore, be presented in a similar way as in the case of the q -independent prescriptions 

BBES, viz., 



transverse phase = V exp 



dr Z^- A^(oo-,0+;M 



(26) 



where I is an arbitrary two-dimensional vector. This result is crucial for our considerations 
in the next section. 



III. CALCULATION OF THE ANOMALOUS DIMENSION IN THE ML-GAUGE 



We start our anomalous-dimension considerations by recalling the modified definition of 
the TMD PDF, proposed in Refs. [li,[l3: 



1 r d^_(f^ ik+^-_^ik±-^± 



It 



2 J 27r(27r)2 

X [oo~, cxD~, oo_l]''^7~'^[oo~, oo_l; oo", 0_l][oo", Oj_; 0~, Oj_] 

x^{0-,Oj_Mp)) R{p+,n-\r,^±) . (27) 

This definition differs from the standard one, given by Eq. ([5]), because it takes into account 
an additional soft factor R which is defined as the vacuum expectation value of the gauge 
links [li, [11 



R 



V exp 



^9 I dC^ eAliO] ■ exp [-^9 t^A'^i^ + Q 



f28) 



illustrated in Fig. [2l and with the involved contours being defined by 

Tcusp : C = {[Pp , -oo < s < 0] U [n's' , < s' < oo] U [Z^r, < r < oo]} 

Tcusp : C/. = , +00 < s < 0] U [n~s' , < s' < oo] U [l^r, < r < oo]} . (29) 

[Note that the result does not depend on the particular choice of the vector l±.] The 



(oo-, 0+, Ix) 



(r. c>o+. 




(0-, -0O+, 0±) 

FIG. 2: The integration contour associated with the additional soft counter term. 

introduction of the soft factor R is necessitated by the demand to cancel undesirable mixed 
rapidity divergences arising in the calculations with light-like quantities [l, [H, 14, 15, [l^. 
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(a) (6) (c) id) 

FIG. 3: Virtual one-loop gluon contributions (curly lines) to the UV-divergences of the TMD PDF 
in the light-cone gauge — graphs (a) and (b). The graphs (c) and (d) are corresponding contributions 
originating from the soft factor R. Double lines denote gauge links. The vertical ones represent 
the transverse gauge links. The Hermitian conjugated diagrams are not shown. 



171 . |18| . Indeed, we have shown in [131, [ij] , using the light-cone gauge with the advanced, 
retarded, or principal-value prescription, that the anomalous dimension entailed by the UV- 
divergences of graph (c) in Fig. [3] (generated by the soft factor R), cancels the corresponding 
contribution of the TMD PDF, given by graph (a) in the same figure. On the other hand, 
pure rapidity divergences still appear in the UV-finite graphs with real-gluon emissions — see 
Fig.H 

In this section, we shall show that the use of the ML-gauge allows one to avoid rapidity 
divergences in the anomalous dimension of the TMD PDF, while preserving at the same 
time the validity (and gauge invariance) of definition fl27j) . Nevertheless, the rapidity di- 
vergences, disentangled from the UV-singularities, are still present in the ML-gauge as well. 
The resummation of them should be pursued by means of the evolution equation which is 
analogous to the CoUins-Soper one. This issue will be considered elsewhere separately. 

Up to the LO in powers of as, the TMD PDF ( 127|) can be cast in the form 

C = /^°^ + f^'' + Oia's) , /(^) = . + £i , (30) 
where we have separated virtual (see Fig. [3]) from real (see Fig. H]) corrections (labeled 
accordingly). The real-gluon terms flU^ do not contain UV divergences and hence will not 
be considered any further. In the tree approximation, one has 

Ci^^k^) = \ I |^e-'=^«-+^fc-^-(p|^(r,^j7^V^(0-,0^)b) 

= 6(1 - x)6^^\k^) . (31) 
The extraction of the UV-singular part of f^^^ proceeds along the lines of our previous 



works, described in [13|, |lj] . In order to isolate the leading-order UV-divergent terms, one 
has to consider the virtual one-gluon contributions depicted in the diagrams of Fig. [31 These 
diagrams amount to 



^(1) 

''vir 

larization with uj = 4 — 2e 



so that = S*^"^ + S'-^^ The quark self-energy diagram (a) gives (in dimensional regu- 



2. [p, as, /i, e) - + 2.^^L - 9 <-F/i J ^2^)a; _ qY(^q2 + ^q) ^^lWJ^2 ^"^^^ 
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/\ 



(a) 



FIG. 4: Real gluon contributions (curly lines) to the TMD PDF in the light-cone gauge using the 
ML pole prescription ("ML gauge"). Double lines denote gauge links. The Hermitian conjugated 
diagrams are not shown. 



with 



(33) 



After some standard calculations, one has for the prescription-independent g^j^^, ("Feynman") 
term 



2e 



(27r)^ (p-g)2(g2 + ^o) p2 
e)r(2-e) 



r(2 - 2e) 



(34) 



The evaluation of the ML-dependent part 



2e 



(27r)^(p-g)2(g2 + ^o) 



ML 



ML 



p2 



(35) 



is more involved. After the transformation of the numerator, one gets 



^ML 



2e 



(27r)^ 



{p - qY 
{p - 



27^ 



ip 



Let us first consider the following integral [2I, [2 
d'^g 1 1 



(27r)-[g+]ML q\p-qf 



da I d(3 



d'^q e 



(g2 + ^0)[g+]MLp2 



(36) 



(271)"^ 

r(6) 1 

(47r)^/2 p+(-p2)e g 



[g+]ML 

rri -2e) 



(37) 



where we have assumed that the direction of the momentum p of the struck quark is purely 
longitudinal, i.e., = 0. At first sight, the above expression seems to have a double pole 
in 1/e. Expanding the F functions, one, however, finds that Eq. fl371) is finite and does not 
contribute any UV singularities. In contrast, the integral with q^ in the numerator is UV 
singular. In order to calculate it, we use the a-representation to obtain 



d'^q g^ 1 _ ^ 

(27r)- [g+]ML q'ip-qf ^ ~2 



^ dx 

X 



dL e^^^P' 



d 



d'^q e 



i{Lq^-2xL{p-q)) 



% J (27r)- [g^ 



ML 



(3J 
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Taking into account that 



d 1 



(39) 



one finds that there are two parts: one proportional to the other to p^. The first part 



vanishes in the final result by virtue of (n 
(137]), gives 

d'^q qf" 1 



{n )^ = 0. Evaluating the second part, using 



(27r)- [g+]ML q'{p - qf (47r)-/2 p+{-p^Y 1(2 - 2e) ' 

This integral contains a single 1/e-pole and thus contributes to the leading UV-singularity. 
In total, we find for diagram (a) 



E(")(p,a,;/i,e) = 4eynman 



Ols_ 

Air 



Cf r(e) -An 



p2 



r(2 - 2e) 



[(1 



4] . (41) 



Extracting the UV divergent terms in the MS-scheme, one gets (after adding the conjugated 
diagrams) : 



An 



1 



1 - 4) - 7e + 47r 



3q;. 



An 



1 



- 7i5 + 47r 



(42) 



This expression makes it apparent that in the ML-gauge the UV-divergent part of the TMD 
PDF (as well as the finite one) do not contain any extra terms of the form Inp"*" which could 
be related to a cusped contour — in contrast to the results obtained using g~-independent 
prescriptions 13|, |l4] • Moreover, one sees that there is no imaginary term, as well, which is, 
however, necessary in order to reproduce the result in covariant gauges. We shall show next 
how this term arises due to the transverse gauge link at light-cone infinity. 



IV. CONTRIBUTION OF THE TRANSVERSE GAUGE LINK AT LIGHT-CONE 
INFINITY 

The path-ordered composite transverse gauge link at light-cone infinity reads 



Pexp 



+ig j dr/^ ■ A-^(cx)~,0+;Z^r + 



Pexp 



-ig I rfrZ^ ■ A^(oo-,0+;/xr) 
'o 



In leading non- vanishing order, the corresponding diagram {h) in Fig. [3] yields 



(43) 



SmlIP./^.^;^) = ~g Cp/i J J 



2n 

q^-l^)+tO {p-qf 



(27r)2 



■ m^{q)A^{q'm 



(44) 



To evaluate this expression, we employ the gluon propagator in the ML-gauge which corre- 
sponds to the correlation function between the longitudinal and the transverse gluon fields: 



(0|A^(g)A^(g')|0) 



q n~ 



(g2 + ^0)[g+]ML 



'-i){2nY5^'\ 



q + q 



(45) 
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Using the explicit form of the ML-gauge field at light-cone infinity (cf. Eq. (125|) ). the trans- 
verse integral can be rewritten in the form 



Taking into account that 

1 



dq+ 



-tq ' oo 



(27r)2 



l^-A^{q)- 



. l- 



(46) 



ML 



V- 



q^q- 



i7i6{q~^q~ 



(47) 



and using the equation (valid in the sense of distributions, see, e.g., Ref. [9| 

-27ri 5(g+) , 



q+ + zO 

one can change variables in the (5-function to obtain 

d'^q 



Sml(P' e) = 



Aq' 



7+(]3 - q) 



(48) 



(49) 



(27r)'^"'^ ' {p-qy g2 ' 
Finally, by taking the sum of the UV-divergent (a) and (6) contributions (Fig. [3]), we find 



.(a+6)UV 



{p,H,as;e) = Cf < - 



1 

4 ~ 2^ 



ITT 

1 

2 



- 7£; 47r 



(50) 



which yields {'-f'^p-f'^ /2p'^ 



^(a+6)UV 



a. 



(p,/i,a^;e) = — Cp 



TT 



1 /3 ZTT , 

7 4-'^|-* + '" 



(51) 



This result resembles what one finds in covariant gauges. After including the mirror contri- 
bution to graph (b) in Fig. [31 one obtains the following expression 



^ML (p,/i,as;e) = — 6] 

TT 
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- 7i^; + 47r 



(52) 



which is analogous to Eq. fl42p and does not contain an imaginary part. Hence, for the 
Mandelstam-Leibbrandt pole prescription, the UV-singular parts of the TMD PFDs repro- 
duce the result obtained in a covariant gauge, where there are no effects from artifacts of 
gauge-contour obstructions (one encounters when using the light-cone gauge in association 
with the advanced, retarded, or principal- value prescription). 



V. EVALUATION OF THE SOFT FACTOR 

To complete our arguments, we have now to verify whether the modified definition (1271) . 



proposed in [13|, [Ij] using a light-cone gauge in conjunction with the advanced, retarded, or 
PV prescription, remains valid in the ML-gauge as well. The main ingredient of this definition 
is a soft factor which was introduced in order to compensate the extra (mixed) UV divergence 
and associated anomalous dimension originating from a cusped contour. However, the latter 



12 



are absent in the ML-gauge, as we have shown above. Therefore, we have to demonstrate that 
the soft factor in this case does not jeopardize Eq. (1271) . In leading order, the UV singularities 
of the soft factor are generated by the self-energy of the light-like gauge link and the one- 
gluon exchanges between the light-like and the transverse gauge link (see diagrams (c) and 
(d) in Fig. [31 respectively). Thus, one has 

<fs'o?t = €U€U0(aD , (53) 

$(0) _ ^ ^(1) _ ^(1) , ^(1) ^(1) _ ^(c) ^{d) 

soft ' ^soft ^soft-virt ~ ^soft-real ' ^soft-virt ^soft-virt ~ ^soft-virt 



rco pa p luj -iq (cr-r) - 

= V'^C'f / da dr -f4 ^ 3 , (54) 

^ Jo Jo J {27rr2q+ q--ql + tO q+ + tOq- ' ^ ^ 

where the vector is chosen to be light-like: = (p"*", 0~, 0^). Due to the relative positions 
of the poles in the Feynman and the ML-denominators, this integral is zero 33|, i-e., both 
poles are on the same side of the g"''-axis: 

= . (55) 

For the same reason, the contribution of diagram (d) in Fig. 2 vanishes as well, entailing 

'^'loft^tirt — 0- ^he other hand, the contribution arising from real gluons, '^'soft-reai; does 
not contain UV-singularities. Hence, ^^^^ reduces to unity, excluding the appearance of 
any contribution to the anomalous dimension of the TMD PDF related to spurious rapidity 
divergences. This result validates Eq. (1271) also for the case of the light-cone gauge with the 
ML-prescription. 

Further, the anomalous dimension of the modified TMD PDF (1271) coincides, therefore, 
with the anomalous dimension of the standard TMD PDF (cf. ([S])) in the light-cone gauge 
with the ML-prescription. Consequently, the renormalization-group properties of the TMD 
PDF are controlled by the following evolution equation 

1 -mod' - ^ f '■^ f dz 



P±{y,q^,a,) = 7ml 5(1 -2/)<5(2)(fc^-g^) + 0(^2), (56) 

where the anomalous dimension coincides with the standard expression, i.e., 

\ d 3 
7ml = --/i— lnSML(as,e) = 7— C'f + O(a^) . (57) 

Z dfj, 4: 11 



VI. CONCLUSIONS 

This work was devoted to the treatment of mixed rapidity divergences of fully gauge- 
invariant TMD PDFs when employing the light-cone gauge in conjunction with the 
Mandelstam-Leibbrandt pole prescription. To this end, we calculated the leading-order 
contributions to the TMD PDF ensuing from virtual gluon corrections. Exactly these terms 
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contain the UV singularities of the TMD PDF and thereby entail its anomalous dimen- 
sion. We have shown by explicit calculation at one loop that, in contrast to other popular 
pole prescriptions, like the advanced, retarded, or the principal-value one, the Mandelstam- 
Leibbrandt prescription possesses the important property that spurious mixed rapidity di- 
vergences, related to obstructions of the gauge contour, are absent. Correspondingly, the 



soft factor, we introduced in [ISj, [IJ] in order to ensure the cancelation of such artifacts to 
the anomalous dimension of the TMD PDF, reduces in this case to unity, thus preserving 
its validity. 

Phenomenologically, the use of the ML pole prescription in the light-cone gauge will facil- 
itate calculations of TMD PDFs in a factorized description of SIDIS cross sections because 
the contributions to the anomalous dimensions from gauge-contour obstructions in the Wil- 
son lines of the TMD PDFs cancel out, making the insertion of a correcting soft factor 
superfluous right from the start. This is also reflected in the evolution behavior of the TMD 
PDFs which is controlled by the standard anomalous dimension one finds in a covariant 
gauge, where anomalous-dimension artifacts are manifestly absent because factorization is 
complete like in collinear factorization. These aspects and the practical analysis of their 
applications will be considered in more detail elsewhere. 
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